Multistability and localization in coupled nonlinear split-ring resonators 
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We study the dynamics of a pair of nonlinear split-ring resonators (a 'metadimer') excited by 
an alternating magnetic field and coupled magnetically. Linear metadimers of this kind have been 
recently used as the elementary components for three-dimensional metamaterials or 'stereometama- 
terials' [N. Liu et al, Nature Photon. 3, 157 (2009)]. We demonstrate that nonlinearity offers more 
possibilities with respect to real-time tunability and a multiplicity of states which can be reached 
by varying the external field. Moreover, we demonstrate almost total localization of the energy in 
one of the resonators in a broad range of parameters. 



Metamaterials are artificially structured composites 
which exhibit electromagnetic properties not available in 
naturally occurring materials. Such structures are largely 
based on subwavelength resonant 'particles' refered to as 
split-ring resonators (SRRs). Scaling down the size of the 
SRRs allows to realize metamaterials up to Terahertz and 
optical frequencies [l|. However, for exploiting the great 
potential of metamaterials for applications, it is desirable 
to change their effective parameters in real time, i.e., to 
achieve real-time tunability. This is also motivated the 
construction of nonlinear SRRs, whose structure is very 
well suited for enhancing nonlinear phenomena [2*]. An 
effective way of constructing an easily controllable non- 
linear SRR is the insertion of a nonlinear electronic com- 
ponent into its slit [1, 0, IE] ■ The arrangement of a large 
number of nonlinear SRRs into a periodic lattice results 
in a nonlinear magnetic metamaterial which is tunable 
by varying the power of the applied field Q . 

It was recently suggested that dimers, comprised of two 
spatially separated SRRs, can be used as elemental units 
for the construction of three dimensional metamaterials 
('stereometamaterials') [7]. The close proximity of the 
SRRs in the dimer results in relatively strong coupling 
between them. A metamaterial comprised of a large num- 
ber of such metadimers can be utilized as a tunable opti- 
cally active medium [8j]. Moreover, if one or both SRRs 
in the metadimer become nonlinear, the metamaterial it- 
self acquires nonlinear properties, essential for real-time 
tunability. It is therefore of great importance to inves- 
tigate the nonlinear properties of the elementary unit of 
such a dimer-based metamaterial, i.e., of the nonlinear 
metadimer, which can be modeled as system of two cou- 
pled nonlinear oscillators. 

We consider an asymmetric metadimer comprised of 
two SRRs which have slightly different slit widths as in 
Fig. 1, dg^i and dg^2, which differentiates their linear 
capacitances Ci and C2, respectively d, [l3|- That in 
turn differentiates the linear resonance frequencies of the 
SRRs through the relation uji ~ l/\/LCi {i — 1,2), while 
the changes in the inductance L and the Ohmic resis- 
tance R of the SRRs are of higher order. The relative 



orientation of the SRRs in the metadimer can be such 
that the SRRs are either narrow-side coupled or broad- 
side coupled (right and left panel of Fig. 1, respectively). 
The nature of the interaction between the SRRs (electric 
or magnetic or both) is determined by the twist angle 
of the slit of one of the SRRs around the H field with 
respect to the other. In both geometries shown in Fig. 
1, that angle is 180° so that the distance between the 
slits is much larger than their widths. Then, the inter- 
action between the SRRs is predominantly magnetic and 
the electric dipole-dipole interactions can be neglected 
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FIG. 1: (color online) Schematic of a broad-side (left panel) 
and narrow-side (right panel) asymmetric metadimer. 

The metadimer is placed in an alternating electromag- 
netic field with the polarization shown in Fig. 1, so that 
only its magnetic component is capable of exciting in- 
duced currents in the SRR rings. In the equivalent circuit 
picture the metadimer is regarded as a pair of periodically 
driven, nonlinear resistor-inductor-capacitor (RLC) oscil- 
lators coupled magnetically through their mutual induc- 
tance M, driven by identical voltage sources. Consider 
a metadimer with the geometry shown in the left panel 
of Fig. 1, for which the magnetic coupling between the 
SRRs is relatively strong. Then, the equations describing 
the dynamics of the (normalized) charge qi , accumulated 
across the slit of the i— th oscillator, reads [l3| 

qi + Ajv/giz + 791 + (dui/dqi) = eg sin(rjr) (1) 
92 + Am9i + 792 + {du2/dq2) = Eq sin(riT), (2) 

where Am = M / L is the magnetic interaction strength. 
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uiQ = ^bj\i02 a characteristic frequency, e = sin(nt) 
is the induced electromotive (emf) force excited in each 
SRR, 7 ~ R\J \JC\CijL is the damping constant, the 
overdots denote derivation with respect to the normalized 
temporal variable r, and 

with 5 = LOxjuji being the resonance frequency mismatch 
(RFM) parameter, which quantifies the asymmetry of the 
metadimer. The average energy, £'(ot, of the metadimer 
can be obtained from the time-average in one period of 
the Hamiltonian 

^ = 0/1 ^ x2 \ (Pi + P2 - '^^MVW-i) + ui + U2, (4) 

where pi = gi + Xufh and P2 = 92 + Aji/gi, with -pi 
and calculated from Eqs. H]) and ([2]). For normaliz- 
ing the earlier equations, we have scaled charge, voltage, 
time, and frequency by Qc, C^c, and wq, respec- 

tively, where Qc = ^G\CiUc, Uc = ^/d^Jd^Ec, with 
Ec a characteristic electric field amplitude. 

The rich dynamical behavior of the nonlinear 
metadimer can be observed in Fig. 2 where a typical 
bifurcation diagram of the currents ii and 12 is shown as 
a function of the driving frequency fl. That diagram can 
be divided into three regions; the region at left, where 
two stable periodic solutions coexist for a wide Q inter- 
val, the region at right where there is a single stable so- 
lution (except in a narrow Q interval), and the chaotic 
region in between separating the previous two ones. Mul- 
tistability exists from = 0.6 to 0.87, where there are 
two stable states with very different energies; the high 
and low energy state with Etot = 1-43 and Etot = 0-014, 
respectively. Importantly, those two states differ consid- 
erably in the distribution of Etot in the two oscillators. 
We define the energy fractions in oscillators 1 and 2 as 
ei = El/ Etot and 62 = E2/Etot, respectively, where Ei 
and E2 are their energies. The high energy state has 
ei — 0.993 and 62 — 0.007 while the low energy state has 
ei = 0.57 and 62 = 0.43. Thus, in the former case, almost 
all the energy is localized in the first of the oscillators. 

The total energy Etot of the two states and the energy 
fractions ei and 62 of the two oscillators as a function 
of n are shown in Figs. 3 and 4, respectively. In that 
specific case, extreme localization occurs in a rather wide 
ri interval, at least from 0.6 to 0.87. If the metadimer 
is initially driven with low frequency it settles into the 
low energy state. As the frequency increases, it passes 
through the point where that state becomes unstable (at 
O = 0.865), and the metadimer suddenly switches to the 
high energy state. We also note that for fl — 0.885 to 
0.895 a stable periodic state coexist with a chaotic state 
(blue in red and green in black). 

In the right-most region of Fig. 2, there is also a nar- 
row O interval (from = 1.02 to 1.05) where multista- 
bility occurs. Outside that interval, i.e., from = 1.05 




FIG. 2; (color online) Bifurcation diagram of the currents ii 
and 12 as a function of the driving frequency Q, for 5 — 0.95, 
7 = 0.01, X = +1/6, A = 0.12, and eo = 0.06. Decreasing 
Q,: ii and 42 correspond to black and red symbols (a and b), 
respectively; increasing fi: ii and 12 correspond to green and 
blue symbols (c and d), respectively. 




FIG. 3: (color online) Total energy of the metadimer Etot 
as a function of the driving frequency Q, for the parameters 
in Fig. 2, and SI decreasing from 1.3 (black symbols - a ); Q 
increasing from 0.6 (red symbols - 6). 



to 1.3 the dimer state has very low energy (see Fig. 
3). However, there are two specific values of ft, i.e., at 
il = 1.08 and 1.025, where Etot is highly localized. At 
those points the energy fractions are ei = 0.02, 62 = 0.98 
and ei = 0.99, 62 = 0.01, respectively. Notice that the 
state at n ~ 1.025 lies in the region of bistability in this 
part of the diagram, which has the lowest energy, while it 
also exhibits localization. The corresponding higher en- 
ergy state at this frequency has ei = 0.46 and 62 = 0.54, 
so that Etot is almost equally distributed in the two os- 
cillators (see also Fig. 4). 

The currents ii and 12 for n = 0.7 (T = 8.98), which 
lies in the region where both multistability and localiza- 
tion occur, are shown in Fig. 5 as a function of t, both for 
the high (Fig. 5a) and the low (Fig. 5b) energy states. 
In the same figures we also plot cos(rir) which is directly 
proportional to the applied magnetic field. The relative 
phase difference between the applied magnetic field and 
the current in each SRR, which is directly proportional 
to its magnetic moment, determines the response of the 
metadimer to that field. We observe that the currents of 
the metadimer in the high energy state have a relative 
phase difference of almost 180° with respect to the ap- 
plied magnetic field, indicating a diamagnetic magnetic 
response. To the contrary, the currents of the metadimer 
in the low energy state are in phase with the applied 
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FIG. 4: (color online) Energy fractions ei and 62 as a func- 
tion of the driving frequency fl for the parameters in Fig. 2. 
Decreasing fl: ei and 62 correspond to black and red symbols 
(a and b), respectively; increasing Q: 6i and 62 correspond to 
green and blue symbols (c and d), respectively. 





FIG. 5: (color online) Temporal evolution of the currents at 
Q = 0.7, for (a) the high energy state; (b) the low energy 
state. The other parameters are as in Fig. 2. 

magnetic field, indicating a paramagnetic response. 

In conclusion, we have revealed a physical mechanism 
of the intrinsic localization of energy in nonlinear mag- 



netic metamaterials in the study of the dynamics of 
coupled sprit-ring resonators. In particular, we have 
found that extreme localization of energy may occur in a 
slightly asymmetric nonlinear metadimer, which also ex- 
hibits multistability of states and chaos for a rather wide 
range of parameters. For a symmetric metadimer (with 
(5 = 1 and other parameters as in Fig. 2), no localized 
states and chaos have been observed, while multistabil- 
ity still occurs. Thus, it seems that a slight asymmetry 
is required for localization to occur. 

We have found that, in the multistability regions, there 
appear two stable states which differ considerably in their 
energies. The magnetic response of the metadimer is ei- 
ther paramagnetic or diamagnetic, depending on the en- 
ergy of the state of the metadimer (low and high, respec- 
tively). The magnetic response of a magnetic metamate- 
rial comprised of such metadimers is determined by aver- 
aging the response of the individual metadimers. When 
the metadimers are in a high energy state, they may 
respond extremely diamagnetically to an applied field. 
Then, the response of the magnetic metamaterial would 
be described macroscopically by a relatively large and 
negative magnetic permeability parameter. 

We believe that this study may be useful for realiz- 
ing strongly nonlinear effects and energy localization in 
metamaterials comprised of a large number of nonlinear 
split-ring resonators where str ongly l ocalized states in the 
form of discrete breathers [H, Hj, El, [iBl , as well as other 
types of localized excitations such as domain walls and 
envelope solitons [13, Ell were predicted theoretically. 
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